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Abstract 

We report a measurement of the first and second moment of the hadronic mass distribution Mx 
in B ^ Xciiy decays for lepton threshold momenta ranging from 0.9 to 1.6 GeV/c in the center 
of mass frame. The measurement uses BB events in which the hadronic decay of one B meson is 
fuhy reconstructed and a charged lepton from the decay of the other B meson is identified. These 
results are obtained from a 140 fb~^ data sample collected near the T(45) resonance with the Belle 
detector at the KEKB asymmetric energy e'^e~ collider. 

PACS numbers: 13.20.He,14.40.Nd 



3 



The heavy quark effective theory (HQET) combined with the operator product expansion 
provides a framework in which inclusive B decay properties can be calculated. In particular, 
the moments of various inclusive quantities in B ^ XcCf decays, such as the lepton mo- 
mentum the 4- momentum transfer squared and the hadronic recoil mass Mx, can be 
related to the non-perturbativeparameters that also appear in calculations of the semilep- 
tonic decay rate (see Ref. [11,0,1^ for calculations of the moments of Mx)- Measurements of 
these moments can thus improve calculations of the semileptonic B decay rate and lead to 
a more precise determination of the magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element Vcb [4]. In this paper, we present a measurement of the first and second 
moment of Mx, (Mx) and (Mx), for lepton momentum thresholds ranging from 0.9 to 
1.6 GeV/c in the center-of-mass (cm.) frame. 

The analysis is based on the data recorded with the Belle detector ^] at the asymmet- 
ric energy e^e~ collider KEKB |6i|, operating at a cm. energy near the T(4S') resonance. 
KEKB consists of a low energy ring (LER) of 3.5 GeV positrons and a high energy ring 
(HER) of 8 GeV electrons. The Belle detector is a large-solid-angle magnetic spectrometer 
consisting of a three-layer silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement 
of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a super-conducting solenoid coil that provides a 
1.5 T magnetic field. The responses of the ECL, CDC (dE/dx) and ACC detectors are 
combined to provide clean electron identification. Muons are identified in the instrumented 
iron flux-return (KLM) located outside of the coil. Charged hadron identification relies on 
the information from the CDC, ACC and TOF sub-detectors. 

The T(45') dataset used for this study corresponds to an integrated luminosity of 140 fb~\ 
or about 152 million BB events. Another 15 fb~^ taken 60 MeV below the resonance are 
used to subtract the non-BB background. Full detector simulation based on GEANT 01 is 
applied to Monte Carlo (MC) simulated events. The size of the MC samples is equivalent 
to about three times the integrated luminosity. The decay B —>■ D*lv is simulated using 
a HQET-based model 0. The ISGW2 model is used for the decays B Diu and 
B — > D**iu. The non-resonant B — »• D^*^Tiii' component is generated according to the 
model of Goity and Roberts [l^ . 

We use BB events in which the hadronic decay of one B meson is fully reconstructed 
(-Stag) and the semileptonic decay of the recoiling B meson (-Bsignai) is inferred from the 
presence of a charged lepton (electron or muon). After selecting hadronic events jlll |. 
we reconstruct i?tag by searching the decay modes B~^ D^*^'^ Tr~^ , D^*^^ , D^*^^af and 
50 _^ D(*)-7r+,L)(*)-p+,L)(*)-a^" ^12]. Pairs of photons satisfying > 50 MeV and 
117 MeV/c^ < 'm{'~f'~f) < 150 MeV/c^ are combined to form 7r° candidates. mesons are re- 
constructed from pairs of oppositely charged tracks with invariant mass within ±30 MeV/c^ 
of the Kg mass and decay vertex displaced from the interaction point. Candidate and 
p° mesons are reconstructed in the tt+tt^ and tt+tt^ decay modes, requiring their invariant 
masses to be within ±150 MeV/c^ of the p mass. Candidate af mesons are obtained by 
combining a p° candidate with a charged pion and requiring an invariant mass between 
1.0 and 1.6 GeV/c^. candidates are searched for in the K~7t~^, K~tt^ti^, K~i:^t:^t:~, 
K^'ir~^TT~ and K^tt^ decay modes. The K~7T^n^ and K^tt~^ modes are used to reconstruct 
mesons. Charmed mesons are selected in a window corresponding to ±3 times the mass 
resolution in the respective decay mode. D*^ mesons are reconstructed by pairing a charmed 
meson with a low momentum pion, D*^ D^tt^ , D'^tt'^. The decay modes D*^ D^tc^ 
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FIG. 1: Reconstruction of the tag-side B meson. The Mbc distribution after requiring AE < 
50 MeV is shown for fully reconstructed (a) (b) decays. The fit uses a Gaussian function 
for the signal and an ARGUS background function 1J| for the combinatorial background. 



and D*^ — D^j are used to search for neutral charmed vector mesons. 

For each Btag candidate, the beam-constrained mass Mbc and AE are calculated [1 



Mbe = V(^wnP-M^ ae = e%-e;,,^, (1) 

where -Ebeam beam energy in the cm. system, and E]^ are the 3-momentum and 

the energy of the Bt^g candidate, respectively. The signal region for i?tag is defined by the 
selections Mbc > 5.27 GeV/c^ and AE < 50 MeV. If multiple candidates are found in a 
single event, the best candidate is chosen based on AE and other variables. Fig. shows 
the Mbc distribution for AE < 50 MeV. The number of the fully reconstructed 
events, obtained by fitting this distribution with the sum of a Gaussian function for the 
signal component and an ARGUS background function [l^], is 76, 155 ±511 (46, 863 ± 374) 
with a purity of the 74% (81%). 

Semileptonic decays of -Bsignai are selected by requiring exactly one identified charged 
lepton among the particles that are not associated to -Stag- In the lepton momentum range 
relevant to this analysis, electrons (muons) are selected with an efficiency of 92% (89%) and 
the probability to misidentify a pion as an electron (a muon) is 0.25% (1.4%) For 



B~^ tags, we require Qe-Qs < 0, where Qe and Qb are the charges of the lepton and of -Btag, 
respectively. 

The charged tracks and neutral clusters which are associated neither with Stag nor with 
the charged lepton are assigned to the hadronic X system. The 4-momenta of the charged 
particles are calculated with the pion mass except for identified kaons. The missing 4- 
momentum in the event is calculated, 

Pmiss = (PLER + PHEr) - PStag - Pi - Px , (2) 

where the indices LER and HER refer to the colliding beams. As only the neutrino in 
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FIG. 2: The Mx distributions for (a) p} > 0.9 GeV/c and (b) p} > 1.6 GeV/c. The data points 
are the T(45) data after subtraction of the non-BB background. The background from secondary 
or misidentified leptons is indicated by the dashed Une. The combinatorial background in the 
-Btag sample is shown as a hatched histogram. 

B Xciiy should be missing in the event, we require the missing mass to be consistent with 
zero, |M2.^J < 2 GeVVc^ 

The neutrino 4- momentum is inferred from the missing momentum, p,^ = {\Puiiss\,Pmiss), 
and the 4-momentum of the X system is recalculated, 

p'x = (Pler + Pher) - Pstag -Pe-Pu ■ (3) 

The Mx resolution (defined as half width at the half maximum) obtained with this procedure 
is about 200 MeV/c^. 

The background consists of the following components: combinatorial background in the 
i?tag sample, background from secondary or misidentified leptons, and background from 
non-BB events. The combinatorial background is dominant. Its shape (in Mx and M]^) is 
determined from the MC simulation, and the Mbc side-band region (5.2 GeV/c^ < Mbc < 
5.27 GeV/c^, \AE\ < 50 MeV) is used to determine its normalization, separately for each 
lepton momentum threshold. The background from secondary or misidentified leptons is 
significant in the low lepton momentum region. Its shape is estimated from the MC simu- 
lation, and the number of well-reconstructed 5tag candidates is used for its normalization. 
The magnitude of the background from misidentified leptons is corrected using the pion 
misidentification rate measured in — > vr+vr" decays in hadronic events. The background 
from non-Si? events is smallest, and it is subtracted using the data taken 60 MeV below the 
T(4S') resonance. Fig. |21 shows the Mx distributions for different minimum lepton momenta 
in the cm. frame, together with the background estimation. 

To correct for the distortions caused by finite detector resolution, the background- 
subtracted Mx and distributions are unfolded using the Singular Value Decomposition 
(SVD) algorithm [17|. The approach consists in describing the detector response by a corre- 
lation matrix between the true and the reconstructed Mx (Mx) values which is determined 
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FIG. 3: The first and second hadronic mass moments, (Mx) and (M^), for different lepton thresh- 
old momenta. The error bars indicate the statistical and total errors. Note that the individual 
moments are highly correlated. All results are preliminary. 



from the MC simulation. This matrix is decomposed into its eigenvalue components, and 
only the contributions leading to a statistically stable solution are kept. For the unfolding 
of the Mx (Mx) distribution, the measured spectrum is divided into 35 (39) 0.1 GeV/c^ 
(0.333 GeV^/c^) wide bins. For the unfolded distribution, 12 (13) bins are used. 

For each minimum lepton momentum, ranging from 0.9 to 1.6 GeV/c, the background- 
subtracted Mx and M|- distributions are unfolded, separately for and tags, and for 
electron and muon events. For each of these four sub-samples, the values of {Mx) and (Mx) 
are obtained by calculating the mean values of the unfolded Mx and Mj- distributions. The 
results obtained by averaging the four sub-sample results are shown in Fig. El and Tables H] 
and|n] All results are preliminary. This procedure has been tested on MC simulated events 
for the full minimum lepton momentum range and no significant bias has been found. 

We consider three sources of systematic error, shown separately in columns three to five of 
Tables E and |TT1 uncertainty related to detector modeling and residual background subtrac- 
tion, uncertainty related to the unfolding procedure, and uncertainty related to the model 
in the MC simulation. The total systematic error quoted on the moments corresponds to 
the quadratic sum of these three components. 

The uncertainty related to the detector modeling and to the background subtraction 
is estimated from the consistency between the four sub-samples, by varying the selections 
used in the analysis (namely the i?tag signal region and the requirement on M^jgg), and by 
varying the background normalization. The uncertainty related to the unfolding procedure 
is estimated by changing the number of eigenvalue components in the unfolding solution. 
The Xc model uncertainty is obtained by varying the fractions of -B ^ D*ih', B —>■ Diu 
and B D**/D^*hiu within ±10%, ±10% and ±30%, respectively, and summing the 
individual variations in quadrature. These ranges of variation roughly correspond to the 
experimental uncertainties 3| in the isospin averaged branching ratios. 

In summary, we have performed a measurement of the first and second moment of the 
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Pmin 


{Mx) 


(GeV/c) 


(GeV/c^) 


0.9 


2.076 ± 0.009 ± 0.010 


1.0 


2.073 ± 0.009 ± 0.010 


1.1 


2.066 ± 0.008 ± 0.008 


1.2 


2.066 ± 0.008 ± 0.006 


1.3 


2.063 ± 0.007 ± 0.006 


1.4 


2.056 ±0.007 ±0.006 


1.5 


2.048 ± 0.008 ± 0.005 


1.6 


2.042 ± 0.008 ± 0.004 



Detector/ Unfolding Xc model 



background 


0.006 


0.008 


0.003 


0.005 


0.008 


0.003 


0.005 


0.005 


0.003 


0.004 


0.004 


0.003 


0.004 


0.005 


0.003 


0.003 


0.004 


0.003 


0.003 


0.003 


0.002 


0.003 


0.001 


0.002 



TABLE I: The first hadronic mass moment {Mx) for different lepton threshold momenta. The 
first error on (Mx) is statistical and the second is total systematic uncertainty, i.e., the quadratic 
sum of the right-most three columns. All results are preliminary. 



Pmin 


(Mi) 


(GeV/c) 


(GeVVc^) 


0.9 


4.334 ± 0.046 ± 0.051 


1.0 


4.312 ± 0.042 ± 0.049 


1.1 


4.285 ± 0.039 ± 0.037 


1.2 


4.258 ± 0.036 ± 0.036 


1.3 


4.247 ± 0.034 ± 0.030 


1.4 


4.214 ± 0.032 ± 0.028 


1.5 


4.166 ±0.033 ±0.035 


1.6 


4.141 ±0.034 ±0.027 



Detector/ Unfolding model 



background 


0.041 


0.027 


0.014 


0.037 


0.030 


0.013 


0.032 


0.017 


0.012 


0.027 


0.021 


0.012 


0.023 


0.016 


0.011 


0.024 


0.009 


0.011 


0.022 


0.025 


0.011 


0.022 


0.010 


0.013 



TABLE IL The second hadronic mass moment (M^) for different lepton threshold momenta. The 
first error on (M^) is statistical and the second is total systematic uncertainty, i.e., the quadratic 
sum of the right-most three columns. All results are preliminary. 



hadronic mass distribution Mx in B ^ Xciiy decays for minimum lepton momenta ranging 
from 0.9 to 1.6 GeV/c in the cm. frame. The results listed in Tables HI and ITTl are compatible 



with recent values from other experiments [20|, |21[. It is expected that this measurement 
(combined with measurements of the semileptonic branching fraction and of the moments 
of the lepton spectrum in B ^ Xciv) will lead to an improved determination of the CKM 
matrix element \Vcb\. 
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